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This paper presents a theoretical investigation of a water- 
augmented turbofan engine, One in which large quantities of sea water 
are injected into the fan duct section. Results indicate that up to 
three or four times dry thrust and propulsive efficiency are obtained 
depending on vessel speed, fan pressure ratio, and amount of water 
injected. Optimum water injection velocity is investigated. Deviations 
from thermal and dynamic equilibrium in the mixing processes are 


investigated with respect to their effect on overall performance. 
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Saturated liquid 
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Gas generator nozzle 
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CHAPTER I 
INTRODUCTION 


With the present day emphasis on nuclear propulsion and its far 
reaching effects on a modern Navy, advances in other types of marine 
propulsion may appear less necessary and certainly less glamorous. 
However, the significance of small to medium displacement craft ina 
limited conflict has been forcibly brought home by the role such craft 
are playing in Vietnam and in other troubled and far-reaching corners 
of the globe. One of the major features of the limited conflict is the 
speed and mobility of participating units, generally with respect to 
land and air forces; however, naval forces must soon follow with faster 
vessels. With existing types of propulsion, speed increases must be 
paid for with dramatic weight increases often outweighing the speed 
increase. Consequently, one of the most recent fields of endeavor in 
marine propulsion has been to develop a system that supplies enough 
power to attain high speeds and yet is lightweight by comparison. This 
paper deals with one of the most recent proposals, the water-~augmented 
turbofan engine. 

Gas turbines are currently being used to drive conventional screw 
propeller vessels and more exotic systems Such as the water-jet driven 
hydrofoil (PCH-2) of the U. S. Navy. However, little serious thought 
has been given to using turbines for jet propulsive power since the high 


exit velocities inherent in gas turbines (even in turbofan engines) 


Le 


produce characteristically low propulsive efficiency and are thus too 
impractical. Conventional thrust augmentation by water injection was 
considered; however, this generally involved injection of water into the 
compressor air flow or into hot turbine exhaust gases. In the case of 
the former, high purity water would be necessary because of the contact 
with vital moving parts. In both cases the high temperatures of the gas 
generator would vaporize most of the water, and the energy required to 
accomplish this would negate the thrust increases (1). A mixed flow 
turbofan where fan duct and gas generator exhaust through the same 
nozzle would have the same deficiency since high temperatures in the 
combined nozzle would again vaporize most of the water. Davison and 
Sadowski (2) proposed a novel method of injecting sea water into the 
relatively cool flow of a fan duct where the fan and gas generator exhaust 
through separate nozzles. As aconsequence, very little water is 
vaporized, and no moving parts are contacted, making the use of sea- 
water feasible. See Figure 1. 

The purpose of this paper is to optimize the design point parameters 
of bypass ratio, water to gas generator ratio, and those velocity, 
temperature, and pressure ratios that directiy affect the performance of 
this system with respect to specific thrust and propulsive efficiency. In 
particular, it is desired to determine to what extent departures from 
thermodynamic and dynamic equilibrium in the mixing processes affect 
performance and what critical limits exist, if any. It would then have to 
be established experimentally that a mixing duct with two-phase flow in 


the actual conditions would satisfy these requirements. 
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This is both necessary and important for two reasons. First, the 
work of Davison and Sadowski (2), which constitutes the only major 
effort on the subject, deals with one particular engine and one fixed 
set of operational parameters. In other words, they have shown it to 
be feasible for one particular engine but have made no effort to find 
under what conditions such a system would be optimum, and such 
information is necessary for the design of a prototype. Second, an 
experimental mixing duct/fan nozzle section should be constructed to 
experimentally verify the theoretical findings of this paper, especially 
with regard to acceptable limits for velocity and temperature lags 
between water and air and also mixing duct lengths necessary to achieve 
them. 

The data obtained in the paper were obtained from a computer 
program designed to analyze the performance of the water-augmented 
turbofan engine while parametrically varying the critical design parameters 
so that graphical and tabular displays of thrust and efficiency variation 
could be obtained and optimum points determined. 

Briefly, the remainder of the thesis will present a description of 
the model used, develop the equations upon which the computer program 
is based, present methods of analyzing the data, and conclude with 
results, discussion, conclusions, and suggestions for further work to 


be done in this area. 
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CHAPTER II 


DESCRIPTION OF MODEL USED 


4,1 Physical: Layout 


Figure 2 is a block-schematic diagram of the Operating system 
depicted in Figure 1. Important points to consider are that the fan duct 
system is completely separate from the gas generator system with 
separate exhausts and that water injection Occurs in the mixing duct 
entrance after the fan, thus contacting nO moving parts. All stations 
are numbered in Figure 2; and when a numerical subscript is involved 
in an equation, the quantity thus subscripted is measured at that station 
in the engine. It should be noted that the figures do not show fuel tanks, 
water injection pumping systems, or any Other cOmponents contained in 


the actual engine. 


2.2 Assumptions 


Because Of the large number of variables in an analysis of this 
type, many simplifying assumptions had to be made to limit computational 
time. The air was assumed to be behaving as a perfect gas with 
constant specific heats. As is common practice in rough design work 
in turbomachinery applications, Cp was assumed constant at one value 
for low temperatures and constant at a second value for high temperatures. 
Gamma was then computed from the appropriate value of Cp. For the 
purposes of this investigation, a plot of Cp versus temperature of the 


gas was made from data in reference (3). The results of this were two 
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values for Cp for the two temperature conditions. For the combustion 
Chamber, which is the transition between the two extremes, an arithmetic 
mean of the two values of Cp was used. This was justified on the basis 
that the plot of Cp versus temperature was roughly a straight line in the 
region of interest. Viscous effects were ignored with the exception of 
the influence of the water on the air in the fan duct. It was assumed 
that no heat was transferred other than in the combustion chamber, and 
no work was transferred other than in the turbine, compressor, and fan. 
Mixing duct area was assumed constant, and it was also assumed that 
the engine would have only converging nozzles. In other words, both 
gas generator and fan exit at atmospheric pressure; however, if this 
produced an exit Mach number greater than one, the values were 
recomputed to exit at a higher pressure with sonic velocity. The last 
major assumption made in the initial analysis was that the air entering 
the system was dry air with no relative humidity effects. For ease in 
calculations one slug of air was assumed to be entering the gas 
generator per unit time. As a consequence, bypass ratio was also the 
mass flow through the fan section per unit time, and an increase in 
bypass ratio keeping other parameters constant was an increase in fan 


section size holding gas generator size constant. 


WA 


CHAPTER Ill 


ANALYSIS OF DRY TURBOFAN 


With zero water injection, -the analysis of the system is developed 


from H-S diagram relationships and basic equations of engineering 


science. This chapter and the equations contained in it are arranged in 


sections similar to those contained in the computer program for the dry 


turbofan engine. 


3.1 Gas Generator Diffuser 


Ambient total pressure is determined from 


yt 
Pro = Pso| L i ——- 


so that defining diffuser efficiency as 


Pr, ~ so 
Fee 


Pro a Pso 


it follows that 


Pr = Peo + Np Pro - Pso| 


Since no work is done in the diffuser, 


lty = 


Pals Tro = ‘Vso | je = Ms | 
and enthalpy is 


= = V7 Tso] a heeca 


where Hprra is the enthalpy of the air at T 
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(1) 


(2) 


(3) 


(4) 


3.2 Compressor 


The ideal (isentropic) temperature at compressor exit for a specified 


Pra 
Pri 





pressure ratio is 





Define compressor efficiency as the ratio of isentropic work input to 


actual work input; thus 


Ij+2 = Ceol Tras - Fr, Wr (7) 
© 


seat 2 | poe H era ' Tso (8) 
Cee 
3.3 Combustor (Burner) 
Recall from Section 2.2 Gag 3 (cee + Cen) 
2.0 
Tr is fixed by maximum allowable turbine inlet temperature. 
Therefore 


Hr: 7 CPs 


Writing the energy equation across the burner, 


nee Tr2 | ha (9) 
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If it is assumed that injection temperature of the fuel is equal to air 
temperature at that point, and the energy needed to vaporize the fuel 


is included by choosing the lower heating value of the fuel, then 
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(10) 
3.4 Turbine 
For the entire system work input must equal work output; therefore, 


H,,, may be calculated from 

















T4 
e « ll 
A, Hr - Hra} = “Me Hre2- He + as is Hre (11) 
and | : . 
ae 413 - Hra (12) 
Vase lrg 7 bb S——— 
Got 


Turbine efficiency is the ratio of actual work output to isentropic work 


output, so that 


Yr Cor 
An expansion whether isentropic or actual is to the same pressure; 
therefore, 
un 
B-\ 


(14) 
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3.5 Gas Generator Nozzle 
Since no work is done in the nozzle, 
Re Gales 


Knowing that the flow exits to ambient pressure fixes the pressure 


ratio Hes) ‘ 
Pra 
Thus 
Y= | 
— | Pss5 e 





{rt 
a | 
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so susye (15) 








Define a nozzle efficiency 
y — Vra - \Vs5 
NM Saree (Set Lo 


Equation (16) may be solved for Tog and exhaust velocity determined 


from 
= a a ne 
VAS =f 2° Cre [Ts - 35 | (17) 
Sound speed is : 
: Qxy = Ve re ee 


Therefore exit Mach number 
— VAS 
= = 18 
Ms ae. (18) 


Note the comments in Section 2.2 concerning procedure when M. =| 


3.6 Fan Diffuser 
Specifying a separate fan diffuser efficiency Ke ep . the equations 
for the fan diffuser correspond to equations (2) through (5) with subscript 


6 substituted for subscript l. 


o./ Pan 


TTT 
TRG 


We , the fan equations are identical to the compressor equations of 


With a corresponding fan pressure ratio 





and fan efficiency, 


Section 3.2 with the subscripts 7 and 6 substituted for the subscripts 


2 and 1 respectively. 


3.8 Fan Mixing Duct 


No work is done in the mixing; therefore, 


on 





Jaen NOZZIE 

The equations of Section 3.5 are valid for the fan nozzle also with 
subscripts 9 and 8 substituted for 5 and 4 respectively and the proper 
fan nozzle efficiency specified. It should be noted at this point that the 


entropy at all stations is calculated from relationships similar to 





_ Tare Prt 
Sark = eo in| i ie R | :. 7 Seess (21) 
where S,,-, is the entropy at Tgog and Ro. 


3.10 Specific Thrust, Propulsive Efficiency, and Specific Fuel 
Con sumption 


Thrust is determined from the momentum flux equation 
ry _ \o 


Specific thrust based on the total mass flow of air entering the system 
is 
TH 
Si = G[me+ me (23) 
Propulsive efficiency is a measure of how well the system is converting 
the change in kinetic energy of the working medium into thrust power of 


the vehicle and can be shown to be 


Zz 


fro 
\" 


+ the[as= Vo] + tre [Van - Vo) (24) 
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Specific fuel consumption is the number of pounds of fuel required to 
produce one pound of thrust continuously for one hour and can be shown 


to be 


SFC = FAR: Mc 7 3600° F (25) 


TR 
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CHAPTER IV 


FIRST WATER INJECTION ANALYSIS 


The injection of water in the mixing duct does not affect the results 
of Sections 3.1 through 3.7. The equations contained in these sections 
are still valid for the case of water injection. The new analysis is 
based on solving the problem in the mixing duct and fan nozzle and 
evaluating the new specific thrust and efficiency from the information. 

In this first water injection analysis it was assumed that none of 
the water in the mixing duct and fan nozzle sections was vaporized. 

The rationale behind this assumption is that, according to reference (2), 


much less than one per cent of the water is vaporized. 


4.1 Definitions 
The following terms will be used throughout the chapters dealing 


with water injection: 


Bypass Ratio (BR) mM e v Wc 


nw / Me 


Water to Gas Generator Ratio (WGR) 


Mw / Me 


Mixture Ratio (MR) 
The mean velocity at a station is defined after Elliot (4) as the 
velocity that the air and water would attain if they were allowed to 


reach equilibrium at that station with no loss in momentum, 


meVa + mw Vw = [eye + MulN (26) 
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thus, 


Mean Velocity (V) = sy) eae MR: Vw (27) 
, + Mi. 


4 .2eaVixincwDwact 


The momentum equation for the assumptions of the paper is 
frie + Aw | AV = -Adp (28) 
The area occupied by the air is WM ¢ 
ss ss 2 /Cr Vo 


The area Occupied by the water is on w/e VV 
WwW VW 





(29) 


Thus, the total flow area is A= wat a... £ MK 
Fi Pave Pw Vw 





Combining equations (29) and (28) and utilizing the fact that the mixing 
duct is assumed to be of constant area, equation (28) can easily be 


integrated and becomes 























— _> |, _MR 
Ve, = Vy T Ps7 <q 1+ mMR]\ Pave Cn Vy (30) 
The energy equation is 
A VJ 2 . " yi 2 
Ne at ee | = MO Ww | w7 7 YW] = 
Zz ra 
(31) 





Equations (30) and (31) have five unknowns: PS8, VA8, VW8, TS8, and 
TW8. With three equations lacking for a solution, it was necessary to 
specificy ratios of unknowns to Serve as equations. These three 


remaining equations were 


Zo 


TSR BWA = TWS / TSA (32) 
VRBWA = VW8/VAS (33) 
PSRS7 > Rane PSH a 


These admirably suited the eventual task of determining the effect 
of deviations from equilibrium in this mixing process on overall 


performance. 


4.3 Fan Nozzle 

The momentum equation cannot now be integrated to yield an equation 
of the form of equation (30) since variation of area and pressure in the 
nozzle is unknown. Several attempts were made to provide a second 
equation either from known facts or from assumptions concerning area 
variation in the nozzle with varying results. The problem was rendered 
academic when the underlying assumption of no water vaporization was 


found to yield erroneous results. (See Section 4.4.) 


4.4 Conclusion of First Attempt 


A continuity equation was introduced into the mixing duct 


calculations in order to reduce the number of variable ratios specified. 


It was introduced in the form of 7 
AN. = Ag 
or more specifically 
Me eee im Me “i Mw (35) 
rs VAq Fwy VW 4 Pr Vas, ears Vwes, 


The problem was now overspecified with more equations than unknowns. 
When two of the three ratios in equations (32)-(34) were used with 
continuity, momentum, and energy as defined in equations (30), (31), 
and (35), ridiculous values were produced for the unknown variables. 
Careful analysis of the equations involved revealed that they were correct 
for the given assumptions indicating that the assumptions needed 
re-evaluation. The net result was that the assumption conceming no 
vaporization of the water was not correct, and, indeed, vaporization Of 
the water was found to play a large role in determining the solution. 
Admittedly a very small percentage of water is vaporized, but for two 
reasons its effect is significant. First, the entalpy of water vapor is 
roughly one-thousand BTU's per lbm. greater than that of liquid water. 
Second, the amount of water present with respect to air (MR) is a number 
ranging in value from roughly ten to four hundred. Thus the contribution 
in an energy equation of vaporized water had to be taken into account 


along with the resultant change in water mass flow. 
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CHAPTER V 
SECOND WATER INJECTION ANALYSIS 


The second analysis is similar to that of Davison and Sadowski (2). 
However, it differs in that their work involves the assumptions of thermal 
and dynamic equilibrium in the mixing duct, dynamic equilibrium in the 
fan nozzle, and zero pressure change in the mixing duct while this 


paper takes into account possible variations in these quantities. 


5.1 Steam Table Properties 

It was necessary in the second analysis to calculate the properties 
of enthalpy of the liquid (Hf), enthalpy change by vaporization (Hfg), 
enthalpy of the saturated vapor (Hg), entropy of the liquid (Sf), entropy 
change by vaporization (Sfg), entropy of the saturated vapor (Sg), and 
partial pressure of the vapor (PPv). These are all complicated functions 
of temperature. A least squares method was used to fit data points 
read in from the steam tables (5) in the range of interest, and it was 
found that a cubic approximation gave accuracy down to the last 
significant figure in the steam tables. The exception was partial pressure 
of the vapor, the actual equation of which was simple enough to be 


included in the program directly. Since 


and 


28 


it was necessary Only to compute the properties Hf, Hfg, Sf, and Sfg, 
thus determining Hg and Sg. The equations generated by the four least 
squares approximations and the PPv calculation are 

HF = -5,2244716x10% + (1.1719422) (Tw) - (3.1882153x10 ~) (tw) 


+(1.9421466x10 ”) (Tw) 3 (38) 


HFG = 1.440548x10° - (1.1017804) (Tw) + (1.0967588x10~°) (Tw) 2 


~(7,4163358x10 ”) (Tw) > (39) 
SF = -1.6947235 + (5.325x10 °) (Tw) - (4.7499752x10 °) (tw)? 
+(1.8842757x10 °) (Tw)> (40) 


SFG = 8.812883 - (2.442480x10 2) (Tw) + (2.8390143x10 °) (Tw) 


~(1, 2441932x107 8) (Tw) ° (41) 


Equations (38) through (41) are valid in the range of temperatures 
from 510 PR to 660 CR. Tw, the temperature of the liquid, is in degrees 
Rankine, and the units of output for the curve fit are BTU/lbm for enthalpy 
and BTU /lbm-°R for entropy. In the equations that follow it will be 


assumed that Hfg, for example, has already been converted into ft-lb/ 


Sug . 
3 
= _ X JA + BX +CX 
where 
Pc = 218.167 atm. 
X= (647.27 - Ty o® 


T = (Tw - 492)(9/5) + 273.16 SK 


ap) 


A = 3.2437814 
7 3 

Pe Gasca rea 

C = 1.1702379x10 © 

D = 2.1878462x10 ° 


and the answer, PPv, will be in atmospheres. Again, in the equations 


that follow, it will be assumed that PPv has been converted into lb/ft? . 


Se aviieing Duct 


The specific humidity, X, of the mixture is defined as the ratio of 
the mass of water vapor to the mass of dry airin the system. From 


perfect gas equations for dry air and water vapor, it can be shown that 


Peve a 
_ {| ——_—_———-_ | /_ Mwy Ise 
X39 ~ Fo, ais Peve MWa | we | (43) 





It was assumed that vapor temperature was equal to liquid temperature. 


Taking into account vaporization, the energy equation is now 


Me [ Msi si Vat 
poe 





+ hay | Hg + Va] = he [Hse Yas] (44) 


o [9 = Xe he] 





2 ; 2 
Hfe + wwe] + Ree | Hae + Vas 


The momentum equation (29) and its integrated result (30) are still 
valid for the case of vaporization. The problem is now one of three 
equations: (44), (43), and (30) with six unknowns: VW8, VA8, TW8, TS8, 
PS8, and X8. By again specifying ratios of unknowns, equations (32) 


through (34), a solution is possible. However, an explicit solution is 


30 


not possible since equations (43) and (44) contain involved functions 
of TW8 through the enthalpy and partial pressure terms. VA8 and VW8 
may be solved using equations (30), (33), and (34). The remaining 
equations may be manipulated in such a manner as to yield a function 
of TS8 only (called FUN) which is equal to zero. A Newton Rhapson 


iterative solution may then be applied 


TS8, = TS8 - FUN/DPUN (45) 


} 


where Bo is the jth approximation, pone is the (j-1)th approximation, 
and DFUN is the derivative of the function with respect to TS8, both 
evaluated at TS8,_y . When the correct value of TS8 is reached, FUN will 
equal zero, and TW8 and X8 may be calculated from equations (32) and 


(43). 


2.3 Fan Nozzle 
Along the same reasoning as in Section 5.2, specific humidity at 


the fan nozzle exit is 


X 2 PPv4 MwWy <a (46) 
4 Fsq al “Pev q MW A Wwv9 


Note again vapor temperature was assumed equal to liquid temperature. 
With varying area, the momentum equation cannot be solved unless 
a geometric configuration is specified. Therefore, the assumption was 


made that the net change in entropy through the nozzle is zero. That is 
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Since exit pressure is atmospheric, the equations (48) and (46) have TS9, 


TW9, and X9 unknown. By specifying 
TSROWA = TW9/TSY (49) 


and solving equations (46) and (48) for a function of TS9Y equal to zero, 
the Newton-Rhapson iterative method again yields values for TS9, TWO, 


and X9. The energy equation may be written 


™M™ = 
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Let VROWA = VW9/VA9 (51) 


and solve equations (50) and (51) for VW9 and VAY 


[Cre 
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2.4 Specific Thrust and Propulsive Efficiency 
Thrust, specific thrust, and propulsive efficiency are defined as 
before, but account must be taken of the varying mass flow rates in the 


mixing sections. 


TH = M+ Vas = Me Vo + [mF + kq the | Vag me Me Vo 


+{thwy - Kame | Vie - ™w1Vwo + As [* ee Po] (5 3) 
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It should be noted that ambient air velocity was assumed equal to ambient 


water velocity. In other words, wind and current were neglected. 


G6 


CHAPTER VI 


RESULTS AND DISCUSSION 


6.1 Organization of Output 

The results obtained focus on two basic indicators of water~augmented 
turbofan performance. These are the ratio of specific thrust for the water- 
augmented engine to specific thrust for the dry engine and the ratio of 
propulsive efficiency for the wet engine to propulsive efficiency for the 
dry engine. Hereafter these two parameters will be referred to as thrust 
ratio and efficiency ratio. The final data concer two broad areas of 
interest, the effect of variations in general system parameters on thrust 
ratio and efficiency ratio holding mixing section equilibrium parameters 
fixed at 1.0, and the effect of variations in mixing section equilibrium 
parameters on thrust ratio and efficiency ratio holding all other system 
parameters constant. 

Subsequent figures show plots of these variations and also the 


effect of water injection velocity on system performance. 


6.2 General System Parameters vs. Thrust and Efficiency Ratios 


The parameters originally thought to have the most significant effects 
on thrust ratio (TR) and efficiency ratio ( 7 R) are water to gas generator 
air ratio (WGR), ship speed (VO), and fan pressure ratio (FPR). Performance 
is also improved by increasing the bypass ratio (BR), but this parameter 
was not thoroughly investigated since its limit is a function of design 


considerations only. 
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Figure 3 presents thrust ratio versus WGR for varying VO and FPR. 
The more water injected, the greater is its influence on the airas 
evidenced by a decrease in fan nozzle exit velocity (VA9). From the 
thrust equation, equation (33), it can be seen that a point is reached 
with increasing WGR where the negative contribution of the momentum 
drag of the water outweighs the positive contributions to thrust. That is, 
the momentum flux into the system eventually increases more rapidly 
than the momentum flux out. Thus, in general, the curves increase, 
reach a maximum or peak value, and then fall off as more water is added. 
Note that for a higher value of VO, values of TR are much lower. This 
is because of the increase in inlet momentum drag of the air and water 
with respect to other thrust contributions. Also note that at the higher 
values of VO, the curves peak at a lower WGR. Again, this is because 
the inlet momentum drag of the water predominates over the increase in 
thrust much more rapidly. 

In comparing Figure 3 with the curves obtained by Davison and 
Sadowski (2), close correlation is found in shape and numerical values 
with differences attributable to the fan pressure ratio. In this paper fan 
pressure ratio is defined as the ratio of total pressure aft of the fan to 
total pressure before the fan, thus making mixing duct static pressure a 
variable with ship speed, while Davison apparently defines fan pressure 
ratio as the ratio of static pressure aft of the fan to ambient pressure, 
since his mixing duct static pressure is constant for all ship speeds. 


Figure 3 shows the fan pressure ratio to be an extremely important design 
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parameter when water injection is considered. The values of thrust 
ratio obtained clearly indicate that excellent values of peak thrust 
augmentation could be obtained in the order of upwards of three times 
dry thrust to over one and one-half times dry thrust. 

Figure 4 presents efficiency ratio versus the same independent 
parameters with a constant fan pressure ratio. The shape of the curves 
obtained is nearly identical to those in Figure 3 with roughly similar 
numerical values. Efficiency ratios peak at slightly higher values of 
WGER than the corresponding thrust ratio indicating that efficiency ratio 
is slightly less sensitive to the increase in momentum drag of the water. 
No particular peak value predominated over the range of velocities for 
either thrust ratio or efficiency ratio. However, for a specified velocity 
and fan pressure ratio, a nearly identical value of optimum WGR was 
found for both thrust and efficiency ratios. Significant was the fact that 
at higher velocities this value was relatively low indicating that high 
speed craft are equally if not more suitable for water injection that those 


operating at lower speeds. 


6.3 Equilibrium Parameters vs. Thrust and Efficiency Ratios 


For the analysis of variations in mixing duct and fan nozzle 
equilibrium parameters on thrust and efficiency ratios, craft speed was 
fixed at fifty knots, fan pressure ratio at 1.5, bypass ratio at 4, WGR 
at 50, and compressor pressure ratio at 13.8. The four equilibrium 
parameters chosen were the ratio of water to air temperature at mixing duct 


exit (TW8/TS8) and fan nozzle exit (TW9/TS9), the ratio of water to air 
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velocity at mixing duct exit (VW8/VA8), and fan nozzle exit (VW9/WAQ). 
Four runs were made holding all but one of the equilibrium parameters 
fixed at 1.0 and varying the remaining One through its maximum possible 
range, (Figures 5 through &). Fora slightly more meaningful output the 
two temperature ratios were varied simultaneously with the velocity 
ratios held constant (Figure 9) and likewise the velocity ratios varied 
together with temperature ratios fixed (Figure 10). 

Figures 5 and 7 indicate that for all other equilibrium parameters 
fixed at unity, equilibrium at the mixing duct exit is either detrimental 
or has no effect on thrust and efficiency ratios. This result was at first 
thought erroneous but then rationalized as being correct (although 
unrealistic). Note that for these curves, equilibrium conditions are 
required at the exit of the fan nozzle. Thus by requiring less equilibrium 
in the mixing duct, the mixing is being forced to take place in the 
nozzle, where it occurs isentropically, resulting in higher efficiency 
and thrust ratios. 

The results obtained in Figures 6 and 8 are much more significant 
as they represent not only physically possible states but, in the case 
of Figure 6, physically probable ones. Figure 6 shows that thermal 
equilibrium at the nozzle exit is definitely not desirable. Peak thrust 
occurs for values of TW9/TS9 between 1.0 and 1.05, which can amount 
to as much as twenty-five or thirty degrees temperature difference. At 
greater values of this temperature ratio, thrust falls off again although 


much more gradually, indicating that a wise operating condition would 
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be to the right of the maximum point. Efficiency ratio is not affected 
to nearly the extent of thrust ratio although it also increases with a 
divergence from equilibrium and in fact never does begin to peak. 
Figure 9 plots thrust and efficiency ratios versugtemperature ratio 
at nozzle exit varying concurrently with temperature ratio at mixing duct 
exit. In this case, the predominance of the nozzle exit temperature 
ratio is again felt as perfect thermal equilibrium is seen to be the worst 
possible condition with a peak value occurring shortly after the 
equilibrium condition. Figure 10 shows the same type of plot with both 
velocity ratios varying and both temperature ratios held constant. It 
reflects the fact that the velocity ratio in the mixing duct makes little 
or no difference on thrust holding the other ratios fixed in that it is 
almost identical to Figure 8. Again departure from equilibrium means 


a decreasing thrust ratio and efficiency ratio. 


6.4 Other Results of Analysis 


One of the more significant results of the paper was discovered 
because Of initial failures to assume realistic data. Initially a constant 
water injection velocity was specified regardless of craft speed. It 
was felt that this could be achieved with proper spray nozzle design and 
a pumping system to provide the necessary pressure differential when 
the speed of the vessel through the water did not create enough head to 
inject the water at the desired velocity. This constant injection velocity 
specification gave absurdly high values of thrust ratio not believed to 


be possible. With the new assumption that the injection velocity is a 
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fraction of the vessel speed (obviating the use of pumps) data were 
obtained which closely correlate with previous work (2). However, the 
fact remained that, fora higher injection velocity, values of specific 
thrust were extremely high. Several TR versus WGR curves for varying 
injection velocity are presented in Figure 11 with equilibrium parameters, 
bypass ratio, and fan pressure ratio fixed. It was then postulated that 
if a pumping system were incorporated to provide the higher injection 
velocities, such a system could be used in emergencies and/or for 
important missions when a burst of additional thrust was required. 
Naturally, the power required to run such a pumping system would reduce 
the propulsive efficiency, but this would be well offset by the otherwise 
unobtainable increase in thrust. Note that this would not require more 
water injection but only higher injection velocities. In effect, it would 


be a type of sea-borne "afterburner." Such a device in intelligence 
ships such as the USS PUEBLO would be of unquestioned value for use 


in emergencies. Perhaps retractable hydrofoils could be incorporated to 


be used only when on the high thrust mode of operation. 
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CHAPTER VII 


CONCLUSIONS AND SUGGESTIONS FOR FURTHER STUDY 


7.1 Conclusions 

The proposed concept of water injection appears to be a feasible 
means of increasing thrust and propulsive efficiency. Peak values of 
thrust ratio and efficiency ratio occur when momentum drag of the 
injected water predominates over the thrust due to additional mass flow. 
These peak values occur at lower values of water to gas generator ratio 
for increasing vessel speed because the momentum drag of the water is 
proportional to vessel speed. 

Peak values for thrust ratios occur slightly before peak values for 
efficiency ratios since the efficiency is less sensitive to inlet momentum 
drag of the water. 

Static pressure in the mixing duct plays an important role in the 
slope of the curves and thus determines the effectiveness of the water 
addition scheme at various ship speeds. 

The equilibrium parameter in the mixing process which appears to 
have significance in affecting system performance is the ratio of water 
temperature to air temperature at the fan nozzle exit. Peak performance 
with regard to thrust ratio is obtained when TW9/TS9 is in the range 
of 1.02 to 1.05. Considerably better values of thrust ratio occur on 
the side away from equilibrium as opposed to closer to equilibrium. 
Efficiency ratio continues to improve as TW9/TS9 moves further away 


from equilibrium. 


47 


Normally, injection velocity would be a fraction of the vessel 
speed determined by how well the ram scoop inlet is designed. Ifa 
pumping system is installed such that the water might be injected ata 
higher velocity, short periods of high thrust could be experienced with 
a corresponding decrease in propulsive efficiency, such as in an 


aircraft afterburner. 


7.2 Suggestions for Further Study 


There is a great deal of further work required before an experimental 
setup can be constructed, much of this work can be accomplished with 
this program and modifications to it. 

The ambient conditions were fixed throughout this analysis and an 
investigation into the variations of system performance for different sea 
water temperatures and air temperatures and pressures could be made. 
Since vaporization of the water in the mixing sections was of such 
significance, the effect of relative humidity in the air should be 
investigated. It should also be possible to choose a friction factor for 
the mixing duct and through the use of the Fanno flow equations 
determine the pressure drop in the mixing duct rather than specify it as 
was done in this analysis. The question of water injection velocity and 
a pumping system "afterburner" should be investigated with regard to the 
decreases in efficiency caused by the additional power required to run 
the pump. 

Last and most important, it remains to design and construct an 


experimental fan section that can investigate the effects of (a) the 
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equilibrium ratios, (b) static pressure in the mixing duct, and (c) injection 
velocities of the air and water on thrust and efficiency and to correlate 
the experimental results with the theoretical values obtained in this 


program and in this analysis. 
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APPENDIX I 


COMPUTER PROGRAM NOMENCIATURE 


Adowill Input Parameters 


Card 1: 


ETAD 
ETAC 
ETAB 
ETAT 
ETAN 
ETAFD 
ETAF 


Card 2: 


DSO 
PSO 
TSO 
HVF 
TREFA 


Card 3: 


AMN7 
TWO 
PRRSZ 


Card 4: 


NVO 
NCPR 
NFPR 
NTT 
NBR 
NWGR 
NMDPR 
NMDVR 
NMDTR 
NFNVR 
NFNTR 


Gas generator diffuser efficiency 
Compressor efficiency 

Bumer efficiency 

Turbine efficiency 

Gas generator nozzle efficiency 
Fan duct diffuser efficiency 

Fan efficiency 


Ambient density (slug/ft?) 

Ambient pressure (lb/ft2) 

Ambient temperature (degrees Rankine) 
Heating value of fuel (BIU/lbm) 
Reference air temperature (492 deg. R) 


Air injection Mach number 
Sea water temperature (deg. R) 
Pressure ratio across burner 


Number of vessel speeds 

Number of compressor pressure ratios 
Number of fan pressure ratios 

Number of turbine inlet temperatures 
Number of bypass ratios 

Number of water to gas generator ratios 
Number of mixing duct pressure ratios 
Number of mixing duct velocity ratios 
Number of mixing duct temperature ratios 
Number of fan nozzle velocity ratios 
Number of fan nozzle temperature ratios 


onl 


Card 5: 


GC 
Costh 


Card 6: 

VO(I) 
Came; 

PTR21 (J) 
Card 8: 

PTR76(M) 
Card 9; 

TT3(K) 
Card 10: 

BR(L) 
Canc: 

WGR(N) 
Cardelz: 

PSR87 (JJ) 
Card 13: 

VR8WA (IK) 
Card 14: 

TSR8WA (IL) 
Cardelo: 

VROWA(IM) 
Card 16: 


TSROWA(IN) 


Cool temperature value of ratio of specific heats 
Hot temperature value of ratio of specific heats 


Values 


Values 


Values 


Values 


Values 


Values 


Values 


Values 


Values 


Values 


Values 


of vessel speed (ft/sec), I= 1,NVO 

of compressor pressure ratio (lb/ft) ; J = 1,NCPR 
of fan pressure ratio (lb/ft) , M=1,NFPR 

of turbine inlet temperature (deg R), K = 1,NTT 

of bypass ratio, L=1,NBR 

of water to gas generator ratio, N = 1,NWGR 

of mixing duct pressure ratio, IJ = 1,NMDPR 

of mixing duct velocity ratio, IK = 1,NMDVR 

of mixing duct temperature ratio, IL=1,NMDTR 
of fan nozzle velocity ratio, IM = 1,NFNVR 


of fan nozzle temperature ratio, IN= 1,NFNTR 


a2 


Al.2 Output parameters 


VA8 


VW8 


VAS 


VW9 


TS8 


TW8 


TS9 


TW9 


STHRUS 


AMNS 


AMN9Y 


ETAPRO 


SFC 


FAR 


XQ 


Air velocity at mixing duct exit (ft/sec) 

Water velocity at mixing duct exit (ft/sec) 

Air velocity at fan nozzle exit (ft/sec) 

Water velocity at fan nozzle exit (ft/sec) 

Air static temperature at mixing duct exit (°R) 
Water static temperature at mixing duct exit (°R) 
Air static temperature at fan nozzle exit (OR) 
Water static temperature at fan nozzle exit (OR) 
Specific thrust (lb thrust/lbm/sec) 

Gas generator exit Mach number 

Fan nozzle exit Mach number 

Propulsive efficiency 

Specific fuel consumption (lb fuel/lb thrust-hr) 
Fuel-air ratio (lb fuel/lb air) 


Specific humidity at fan nozzle exit (lb vapor/lb dry air) 
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